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Abstract 

We analyse the effects of thermal conduction in a relativistic fluid, 
just after its departure from hydrostatic equilibrium, on a time scale 
of the order of thermal relaxation time. It is obtained that the result- 
ing evolution will critically depend on a parameter defined in terms of 
thermodynamic variables, which is constrained by causality require- 
ments 
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1 Introduction. 



Most part of the life of stars (at any stage of evolution), may be described on 
the basis of the quasi-static approximation (slowly evolving regime). This is 
so, because most relevant processes in star interiors take place on time scales 
that are usually, much larger than the hydrostatic time scale 
However, during their evolution, self-gravitating objects may pass through 
phases of intense dynamical activity for which the quasi-static approxima- 
tion is clearly not rehable (e.g., the quick collapse phase preceding neutron 
star formation). All these phases of star evolution ("slow" and "quick") are 
generally accompanied by intense dissipative processes, usually described in 
the diffusion approximation. This assumption, in its turn, is justified by 
the fact that frequently, the mean free path of particles responsible for the 
propagation of energy in stellar interiors is very small as compared with the 
typical length of the star. 

In this work we shall study the influence of thermal conduction on the evolu- 
tion of a self-gravitating system out of hydrostatic equilibrium (in the "quick" 
phase). 

However, instead of following its evolution long time after its departure from 
equilibrium, we shall evaluate the system immediately after such departure. 
Here "immediately" means on a time scale of the order of thermal relaxation 
time, before the establishment of the steady state resistive flow. 
Doing so we shall avoid the introduction of numerical procedures which might 
lead to model dependent conclusions. 

On the other hand, however, we shall obtain only indications about the 
tendency of the object and not a complete description of its evolution. 
As we shall see, there appears a local parameter formed by a specific combina- 
tion of thermal relaxation time, thermal conductivity, proper energy density 
and pressure, which critically affects the evolution of the object and which 
is constrained by causality requirements. 
The paper is organized as follows. 

In the next section the field equations, the conventions and other useful 
formulae are introduced. In section 3 we briefly present the equation for 
the heat conduction. The central problem is analysed in section 4 and a 
discussion of results is given in the last section. 



2 



2 Field Equations and Conventions. 

We consider spherically symmetric distributions of collapsing fluid, which for 
sake of completeness we assume to be anisotropic, undergoing dissipation in 
the form of heat flow, bounded by a spherical surface S. 
The line element is given in Schwarzschild-like coordinates by 

ds^ = e^df - e^dr^ - (dO'^ + sin^ed(p'^) (1) 

where v{t^r) and A(t,r) are functions of their arguments. We number the 

coordinates: x° = t; = r; = 6; = 0. 

The metric ([|) has to satisfy Einstein fleld equations 

= -SttT; (2) 

which in our case read f^: 

-8-T/ = -l + e-(l + ^) (4) 



- SvrTi = -SttTI = - ^ (2A + A(A - u] 

+ ^ (2." + - AV + (5) 

-87rToi = -- (6) 
r 

where dots and primes stand for partial differentiation with respect to t and 
r respectively. 

In order to give physical signiflcance to the TjJ' components we apply the 
Bondi approach 0]. 

Thus, following Bondi, let us introduce purely locally Minkowski coordinates 
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dr = e'^/^dt dx = e^/^dr 



dy — rdO dz — rsin6d(f) 



Then, denoting the Minkowski components of the energy tensor by a bar, we 
have 



-t n 



T° Tl = Tl r| = T| T| = T| Toi = e'^^'+^^'^T^ 
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Next, we suppose that when viewed by an observer moving relative to these 
coordinates with velocity uo in the radial direction, the physical content of 
space consists of an anisotropic fluid of energy density p, radial pressure P^, 
tangential pressure and radial heat flux q. Thus, when viewed by this 
moving observer the covariant tensor in Minkowski coordinates is 



/ p -g \ 

-q Pr Q 

Pi 

V Px / 



Then a Lorentz transformation readily shows that 



rp\J rp\) 



puj 



(l-a;2)i/2 



Tn2 

-^2 



m3 



rp2 

^2 



(7) 

(8) 
(9) 



(p + Pr)a;e 



(i^+A)/2 



;//2 A 



(l-u;2-)i/2 



[l + uj') (10) 



with 



ge 



-A/2 



1-6^2)1/2 



(11) 
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Note that the velocity in the {t, r, 6, 0) system, dr/dt, is related to u by 

dt ^ ' 

At the outside of the fluid distribution, the spacetime is that of Vaidya, given 
by 

^{^_ j + 2dudR - (de'' + sm^Odct)^) (13) 

where m is a time-like coordinate such that u = constant is (asymptotically) 
a null cone open to the future and i? is a null coordinate {gnB. = 0). It should 
be remarked, however, that strictly speaking, the radiation can be considered 
in radial free streaming only at radial infinity. 

The two coordinate systems {t, r, 9, (p) and (u, R, 9, 0) are related at the 
boundary surface and outside it by 

u = t-r-2Mln(—-l) (14) 
V2M / ^ ^ 

R = r (15) 

In order to match smoothly the two metrics above on the boundary surface 
f = ^T.{t), we have to require the continuity of the first fundamental form 
across that surface. As result of this matching we obtain 

[Pr]^=\Qe'/^l-u'Y^'] =[g], (16) 



s 

expressing the discontinuity of the radial pressure in the presence of heat 
flow, which is a well known result [Q. 

Next, it will be useful to calculate the radial components of the conservation 
law 

= (17) 

After tedious but simple calculations we get 

(-80' = {Ti - Tf) + 4,.' (r; - r„") + (a + f - f] (18) 
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which in the static case becomes 

representing the generahzation of the Tolman-Oppenheimer-Volkof equation 
for anisotropic fluids [Q. 



3 Heat Conduction Equation. 



As we mentioned in the introduction, in the study of star interiors it is 
usually assumed that the energy flux of radiation (and thermal conduction) is 
proportional to the gradient of temperature (Maxwell-Fourier law or Eckart- 
Landau in general relativity). 

However it is well known that the Maxwell-Fourier law for the radiation flux 
leads to a parabolic equation (diffusion equation) which predicts propagation 
of perturbation with infinite speed (see [^]-|^ and references therein). This 
simple fact is at the origin of the pathologies found in the approaches of 
Eckart [|lOl and Landau fill for relativistic dissipative processes. 
To overcome such difficulties, different relativistic theories with non- vanishing 
relaxation times have been proposed in the past |]l2l-[jl5l. The important 



point is that all these theories provide a heat transport equation which is 
not of Maxwell- Fourier type but of Cattaneo type 0], leading thereby to a 



hyperbolic equation for the propagation of thermal perturbation. 
Accordingly we shall describe the heat transport by means of a relativistic 
Israel-Stewart equation |§] , which reads 



r^ + g" = /tP"^(T^-Ta^)-™°g^a^-l/€T2f^M^) (20) 

JJS 2 / ;/3 

where k, t, T, and denote thermal conductivity, thermal relaxation 
time, temperature, the heat flow vector and the components of the four 
acceleration, respectively. Also, P"^ is the projector onto the hypersurface 
orthogonal to the four velocity m". 

In our case this equation has only two independent components, which read, 
for a = 



6 



re 



{Quj + Quj + QuXj + T Q'u? + QujJ + 



1/2 



KTue-^/^ 1 



2(1 -(^2)1/2 2 



+-rQuj 



UOLO A 



+ rg 7^ 1 X 



(21) 



and for q; = 1 



re 



+ geV2 (l-c.^) = 



(l-a;2)i/2 (l-a;2) 



2^1/2 



+lrQ 
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ujuo' 



where the expressions 



= g (cue(^-^)/2^ 1, 0, O) (24) 

have been used. 

We are now ready to get into the central problem of this work. 



4 Thermal Conduction and Departure from 
Hydrostatic Equilibrium. 

Let us now consider a spherically symmetric fluid distribution which initially 
may be in either hydrostatic and thermal equilibrium (i.e. u; = Q = 0), or 
slowly evolving and dissipating energy through a radial heat flow vector. 
Before proceeding further with the treatment of our problem, let us clearly 
specify the meaning of "slowly evolving". That means that our sphere 
changes on a time scale which is very large as compared to the typical time 
in which it reacts on a slight perturbation of hydrostatic equilibrium. This 
typical time is called hydrostatic time scale. Thus a slowly evolving system 
is always in hydrostatic equilibrium (very close to), and its evolution may be 
regarded as a sequence of static models linked by (^). 

As we mentioned before, this assumption is very sensible, since the hydro- 
static time scale is usually very small. 

Thus, it is of the order of 27 minutes for the sun, 4.5 seconds for a white 
dwarf and 10~^ seconds for a neutron star of one solar mass and 10 Km 
radius 0. 

In terms of u and metric functions, slow evolution means that the radial 
velocity u measured by the Minkowski observer, as well as time derivatives 
are so small that their products and second order time derivatives may be 



8 



neglected (an invariant characterization of slow evolution may be found in 
Thus 10 

j) ^ \ ^ \i/ ^ X"^ ^ i)"^ ^ cu"^ ^ iu = (25) 

As it follows from (|^) and ([T0|), Q is of the order 0{uj). Thus in the slowly 
evolving regime, relaxation terms may be neglected and (^) becomes the 
usual Landau-Eckart transport equation |T^. 

Then, using (^) and (|l^) we obtain ([l9|), which as mentioned before is the 
equation of hydrostatic equilibrium for an anisotropic fluid. This is in agree- 
ment with what was mentioned above, in the sense that a slowly evolving 
system is in hydrostatic equilibrium. 

Let us now return to our problem. Before perturbation, the two possible 
initial states of our system are characterized by: 

1. Static 

dj = Q = uj = Q = (26) 

2. Slowly evolving 

u = Q = (27) 
Q ^ 0{uj) 7^ (small) (28) 

where the meaning of "small" is given by (pSj). 

Let us now assume that our system is submitted to perturbations which force 
it to depart from hydrostatic equilibrium but keeping the spherical symmetry. 
We shall study the perturbed system on a time scale which is small as com- 
pared to the thermal adjustment time. 

Then, immediately after perturbation ("immediately" understood in the 
sense above), we have for the first initial condition (static) 

uj = Q = (29) 

u^Q^O (small) (30) 
whereas for the second initial condition (slowly evolving) 

Q ^ 0(u) ^ (small) (31) 
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Q^u^O (small) (32) 

As we shall see below, both initial conditions lead to the same final equations. 
Let us now write explicitly eq.(|TBp. With the help of (0)-(|TUp, we find after 
long but trivial calculations 



+ T^- ^ + ^ + 



l-UJ^ l-UJ^ l-UJ^ {l-uj'^f 
^ (1 - UJ^f ^ (1 - U^)'/' ^ (1 - ^ (1 - ^2)3/2 



+ — 77^ \ (Pr+pUJ^](p+Pr) n 

+ — (P, + pcj=^ — ^ ^ + — p + P, ^ 



A + ^-^) (33) 



87rr V 2 2 ^ 
which, when evaluated immediately after perturbation, reduces to 



{P + Pr)m ^ AtTT . ^ _^ 2 (P, - P^) ^ 6^ ■ 

r2(l — 2m/r) (1 — 2m/r) V 7- 87J-7- 



for both initial states. 

On the other hand, an expression for A may be obtained by taking the time 
derivative of (IBD 



A = -87rre(^+^)/2[(p + P.)-^!j + geV2. ^ + ^ 



l-a;2 2 ^ n -,^2^1/2 2 
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(p + Pr)u; X 

+ (p + Pr) dj 



to 



which, in its turn, when evaluated after perturbation, reads 

replacing A by ( pUf ) en (p^), we obtain 

- e("-^)/2^ = (p + P,) cI; + 
where i? denotes the left-hand side of the TOV equation, i.e. 



2.2 + 



(1 -Cj2)l/2 



(35) 



(36) 



(37) 



^ _ dPr AirrP^ Priri AnrpPr 

dr 1 — 2m/r r'^{l — 2m/r) 1 — 2m/r 

^ pm 2 (Px - Pr) 

(1 — 2m/r) r 

= P; + ^(p + P.)--(i^±-i^r) (38) 

2 r 

The physical meaning of R is clearly inferred from (^). It represents the 
total force (gravitational + pressure gradient + anisotropic term) acting on 
a given fluid element. Obviously, P > 0/P < means that the total force is 
directed inward/ outward of the sphere. 

Let us now turn back to thermal conduction equation (pO]). Evaluating its 
t-component (given by Eg . (pT|) ) immediately after perturbation, we obtain 
for the first initial configuration (static), an identity. Whereas the second 
case (slowly evolving) leads to 

u (t' + j = (39) 

which is to be expected, since before perturbation, in the slowly evolving 
regime, we have according to Eckart-Landau (valid in this regime) 
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Q = -f^e-' (t' + (40) 

Therefore, the quantity in bracket is of order Q. Then immediately after 
perturbation this quantity is still of order 0{uj), which implies (|39|) . 
The corresponding evolution of the r-component of the equation (^) yields, 
for the initially static configuration 

= -KTdj (41) 
where the fact has been used that after perturbation 

Q = =^ T' = -^ (42) 

For the second case, the r-component of heat transport equation yields also 
(^rp, since after perturbation the value of Q is still given by (^), up to 0{uj) 
terms. 

Finally, combining (|37|) and (p]) 

e('^-^)/2i? 1 ^ ^ 

^ = --, ;tv X -/ — r (43) 

or, defining the parameter a by 

a = —r^ ^ (44) 

-e^-^^/^R=(p + p^)dj{l-a) (45) 

Let us first consider the a = case. Then, last expression has the obvious 
"Newtonian" form 

Force = mass x acceleration 
since, as it is well known, [p + P^) represents the inertial mass density and 
by "acceleration" we mean the time derivative of cu and not (a^a^)^/^. In this 
case (a = 0), an outward /inward acceleration (cj > 0/cj < 0) is associated 
with an outwardly / inwardly {R < 0/R > 0) directed total force (as one 
expects!). 

However, in the general case {a ^ 0) the situation becomes quite different. 
Indeed, the inertial mass term is now multiplied by (1 — a), so that if a = 1, 
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we obtain that u; 7^ even though R = 0. Still worse, if a > 1, then 
an outward /inward acceleration is associated with an inwardly /outwardly 
directed total force!. However as we shall see in next section, causality require- 
ments constrain a to be less than 1. 

The last term in (^) is frequently omitted (the so-called "truncated" the- 
ory) |jT9[. In the context of this work both components of this term vanish 
and therefore all results found above are independent of the adopted theory 
(Israel-Stewart or truncated). 



5 Discussion. 



Restrictions based on stability and causality were derived in ^ for Israel- 
Stewart thermodynamics. According to equation (134) of ||^, it follows that 
we must have a < 1 in order to guarantee that thermal pulses are subluminal 
In fact there is a similar parameter in the case of bulk viscous perturbations 
[^], which due to causality an stability limits should be smaller than one. 
Before concluding we would like to make the following remarks: 

1. Observe the formal similarity between the critical point and the equa- 
tion of state for an inflationary scenario (p = — P, ). 

2. It should be clear that in the context of the perturbation scheme used 
here, we get information only about the tendency of the system. To 
find out the real influence of the critical point on the evolution of the 
object, the full integration of the equations is required. Calculations 
involving such integrations have been performed in the past [^ -[p^ 



however in neither one of the examples examined there, the system 
reaches the critical point. Furthermore, our configurations are initially 
in global thermal equilibrium, which is a highly idealized situation. In 
this sense, it should be stressed that our aim here is not to model a 
real star but to study some specific aspects of relativistic diffusion. 

It should be clear that the analysis presented here depends strictly on 
the validity of the diffusion approximation, which in turn depends on 
the assumption of local thermodynamical equilibrium (LTE). There- 
fore, only small deviations from LTE can be considered in the context 
of this work. 
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4. For the sake of completeness we have considered an anisotropic fluid 
(instead of an isotropic one), leaving the origin of such anisotropy 
completely unspecified. As it is apparent, anisotropy does not affect 



the most important result obtained here (Eq.(^5D). However, should 
anisotropy be related to viscosity, then for consistency the anisotropic 
pressure tensor should be subjected to the Israel-Stewart causal evolu- 
tion equation for shear viscosity 
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